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Abstract:  The  influence  of  an  S  adlayer  on  CO  adsorption  onto  Ni(100)  is  examined. 
Tight  binding  extended  Huckel  calculations  on  a  three  layer  model  slab  indicate  that 
the  interadsorbate  separation  distance  determines  not  only  the  mechanism,  but  also  the 
effect  of  the  interaction.  If  the  C-S  distance  is  short,  sulfur  induces  site  blockage  of  CO 
chemisorption  by  means  of  a  direct,  repulsive  interadsorbate  mechanism.  If  the  separation 
is  increased  beyond  the  normal  S-C  bond  range,  the  sulfur  adatoms  work  indirectly  via 
modification  of  the  electronic  structure  of  the  substrate.  This  is  a  form  of  through-bond 
coupling.  It  is  consistent  with  the  well-documented  sulfur  poisoning  of  CO  adsorption 
and  its  usual  explanation  via  relative  electronegativities  of  adsorbates,  but  there  are  some 
conceptual  differences.  At  longer  coadsorbate  separations,  there  is  an  interesting  reversal 


of  the  bonding  trends,  which  has  some  experimental  support 
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The  effect  of  atomic  adsorbates  on  the  chemisorption  of  small  molecules  and  on  the 
rate  of  certain  catalytic  reactions  is  dramatic1.  Electronegative  impurities,  such  as  halo¬ 
gens  and  chalcogens,  tend  to  hinder  these  processes.  In  contrast,  alkali  metals  and  other 
electropositive  elements  can  function  as  promoters.  A  great  variety  of  experimental  and 
theoretical  studies  have  been  conducted  with  the  motivation  of  understanding  the  poison¬ 
ing  and  enhancement  effects  of  coadsorbates.  For  example,  the  adsorption  of  CO  onto 
sulfided  surfaces  is  often  used  as  a  model  for  the  sulfur  poisoning  of  Fisher-Tropsch  hydro¬ 
carbon  catalysis  from  CO  and  Hj.  We  will  focus  on  S/CO  coadsorbates  on  nickel  surfaces, 
as  the  body  of  experimental  work  dedicated  to  these  systems  is  voluminous. 

As  is  true  of  other  coadsorbate  systems,  the  basic  nature  of  the  interactions  between 
surface  species  is  a  matter  of  controversy.  Opinions  are  substantially  polarized  between 
two  extremes.  The  mechanism  is  described  either  as  dominated  by  delocalized  long-range 
effects  which  allow  a  single  impurity  adatom  to  modify  many  adsorption  sites*,  or  alter¬ 
natively  as  mediated  by  local  bonding  and  short-range  site  blockage5.  Using  the  army 
of  acronymical  methodologies  available  to  surface  science,  experimental  evidence  can  be 
found  to  support  either  side.  We  list  only  a  few  examples.  Goodman  and  Kiskinova2,4, 
and  Erley  and  Wagner5  advocate  the  long-range  theory  based  on  the  well-documented 
nonlinearity  of  both  the  CO  saturation  coverage  and  CO/Hj  catalytic  methanation  as  a 
function  of  preadsorbed  sulfur  on  Ni(100)  and  Ni(lll).  On  Ni(100),  a  one  fifth  sulfur 
monolayer  causes  am  order  of  magnitude  decreases  in  the  rate  of  methanation.  Goodman6 
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concludes  that  each  sulfur  atom  affects  some  ten  Ni  surface  atoms.  Madix  et  al.3,7,8  im¬ 


plicate  local  site  blockage  to  explain  similar  results.  Gland  et  al.9  favor  the  short-range 


model  for  Ni(100)  in  light  of  high  resolution  electon  energy  loss  spectroscopy  (HREELS) 


and  temperature  programed  desorption  (TPD)  studies.  The  infrared  reflection-adsorption 


spectroscopy  (IRAS)  work  of  Trenary  et  al.10  on  Ni(lll)  is  in  good  agreement  with  the 


latter. 


Theoretical  treatments  of  adatom  poisoning  and  promotion  of  CO  chemisorption  are 


generally  presented  in  the  framework  of  the  Blyholder  model11.  The  adsorption  geometry 


is  widely  accepted  to  be  through  the  carbon  end,  exactly  or  nearly  perpendicular  to  either 


the  clean  or  preadsorbed  nickel  surface17*''1.  Regardless  of  the  specific  adsorption  site,  the 


chemisorptive  bond  in  the  Blyholder  model  results  from  electron  donation  from  the  CO 


5a,  la,  into  the  empty  surface  levels  and  backdonation  from  the  surface  into  the  CO  2ir *, 


lb.  As  both  CO  levels  are  localized  on  the  carbon  end13,  they  are  at  least  geometrically 


well  suited  for  the  electron  transfers.  Sulfur  and  other  electronegative  adatoms  withdraw 


\f 


electron  density  from  the  surface.  Little  change  in  5 a  to  metal  donation,  2a,  is  presumed, 
but  metal  to  2n *  backbonding,  2b,  should  be  reduced,  and  hence  the  poisoning  effect.  As 
the  occupation  of  the  C-0  antibonding  2 r*  decreases,  the  C-0  bond  strengthens  at  the 
expense  of  the  chemisorptive  bond.  Electropositive  adatoms  should  work  in  the  opposite 
maimer.  Electron  density  is  pushed  into  the  surface,  thus  enhancing  its  backbonding 
capacity. 

Theoretical  studies  are  able  to  reproduce  poisoning  and  promotion,  and  most  indicate, 
at  minimum,  a  rudimentary  agreement  with  this  simple  MO  picture.  A  cornucopia  of 
calculation  methodologies  has  been  called  upon,  including  LCAO-type  on  clusters140-**  and 
monolayers15,  monolayer  SLAPW-type18,  cluster  LCGTO-Xcr170*,  the  effective-medium 
theory  applied  to  jellium  surfaces1*0’*,  the  muffin-tin  approximation  applied  to  clusters19 
and  moat  recently,  all-electron  local-density-functional  theory20.  Nevertheless,  there  is  no 
more  agreement  (as  to  the  range  or  detailed  nature  of  the  interactions)  between  the  various 
calculations  than  that  found  in  experimental  studies.  Cluster  calculations  are  inherently 
sensitive  to  short-range  effects  and  may  have  difficulty  modeling  the  metallic  nature  of  the 
surfaces.  However,  even  the  conclusions  drawn  from  the  extended  structure  calculations 
alone  are  divergent.  Benzinger  and  Madix15  suggest  that  direct  interadsorbate  interactions 
may  be  the  key  factor,  although  the  metal-mediated  electron  transfers  predicted  from  the 
Blyholder  mechanism  are  clearly  reproduced.  On  the  other  hand,  N0rskov  et  al.18°  *  find 
long-range  electrostatic  forces  to  be  more  significant.  The  FLAPW  results  of  Wimmer  et 


al.20  straddle  the  contoversy.  A  long-range  interaction  is  dominant  for  a  K  adatom,  but 
for  S,  both  direct  and  substrate  mediated  effects  are  surmised. 

In  this  contribution,  we  examine  the  S  +  CO  system  on  Ni(lOO)  with  the  extended 
Huckel  tight  binding  method214’*.  The  procedure  has  well-known  deficiencies  —  for  in¬ 
stance,  it  does  not  give  reliable  potential  energy  curves.  But  it  is  a  transparent  methodol¬ 
ogy  and  reveals  clearly  the  basic  interactions  that  are  responsible  for  bonding. 

We  will  assemble  the  Ni/S/CO  system  from  its  various  component  parts.  Computa¬ 
tional  details  are  given  in  the  Appendix.  As  the  focus  or  previous  work  by  this  group  has 
been  the  adsorption  of  CO  onto  Ni(100)  and  other  surfaces22,  several  of  the  closely  related 
systems  —  the  bare  Ni(100)  surface,  the  unadsorbed  CO  net,  the  Ni(100)-c(2x2)CO  on-top 
and  Ni(lll)CO  bridging  —  are  well  characterized  and  our  discussion  will  be  appropriately 
brief. 

For  the  Ni(100)  substrate,  a  three  layer  slab  with  the  bulk  nickel  lattice  constant 
(a=3.524A23)  was  used.  The  nearest  neighbor  contact  is  2.49A  on  this  square  net  surface. 
The  calculated  density  of  states  (DOS)  exhibits  a  compact  d  band  between  —12  and 
— 8eV,  and  very  disperse  s  and  p  bands  between  —12  and  SeV22,24.  As  the  surface  atoms 
have  a  lower  coordination  number,  these  states  are  less  disperse  than  the  bulk  states. 
Consequently,  if  the  Fermi  level  falls  above  the  midpoint  of  the  d  block,  as  for  nickel,  the 
surface  is  negatively  charged  with  respect  to  the  bulk. 

Whether  the  three  layer  slab  is  an  appropriate  model  for  the  surface  can  estimated 


from  the  charge  on  the  middle  layer.  An  infinitely  deep  system  will  have  a  bulk  charge 
approaching  zero.  With  a  55  k  point  set38,  we  compute  an  excess  charge  of  +0.180e~  per 
atom  for  three  layers,  +0.166e~  for  four  layers  and  +0.157e"  (middle  layer)  for  five.  As  the 
convergence  is  slow  in  this  range  of  reasonably  sized  unit  cells,  we  make  our  choice  on  the 
basis  of  computational  economics.  The  three  layer  substrate  was  used  for  all  calculations. 

Our  first  priority  in  the  study  of  the  S/CO  system  on  Ni(lOO)  is  to  establish  proper 
adsorption  geometries.  So  let  us  review  briefly  what  as  known  from  experiment.  The 
preference  for  atomic  adsorption  at  the  highest  coordination  site  is  nearly  universal  on 
clean  or  coabsorbed  transition  metal  surfaces26.  Specifically,  chalcogen  adsorption  at  the 
Ni(100)  four-fold  hollow  is  well  documented270-*.  CO  chemisorbs  molecularly  and  does 
not  dissociate  to  any  measurable  degree  on  late  transition  metal  surfaces28.  On  clean 
Ni(100),  CO  site  preference  is  coverage  and  temperature  dependent29®-*.  At  saturation 
coverage  and  below  150K,  a  compression  structure  of  both  on-top  and  bridging  CO  is 
formed.  When  the  compression  structure  is  relaxed  to  a  9= j,  c(2x2)  monolayer  (or  if 
coverage  is  further  reduced),  only  the  on- top  mode  is  obseved.  The  HREELS  vibrational 
spectrum9  is  characterized  by  a  single  C-0  stretching  frequency  at  2005cm-1.  If  sulfur  is 
added,  and  its  surface  concentration  is  increased,  this  dominant  peak  is  downshifted  and 
two  new  CO  modes  grow  in,  one  at  a  lower  frequency  and  the  other  at  a  higher  frequency. 
The  adsorption  geometry  of  the  new  modes  is  unknown.  Using  partial  thermal  desorption, 
the  high  frequency  peak  is  correlated  to  low  temperature  desorption,  the  low  frequency 
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to  high  temperature  and  the  dominant  mode  to  an  intermediate  temperature.  The  CO 
layer  is  completely  desorbed  below  430K,  which  is  the  predominant  desorption  temperature 
from  the  clean  surface.  The  inverse  relationship  between  the  C-0  stretch  frequency  and 
desorption  temperature  is  consistent  with  the  greater  susceptibility  of  the  chemisorptive 
bond  strength  to  the  ir,  rather  than  the  a  system,  as  advocated  in  the  Blyholder  mechanism. 
However,  we  expect  a  unidirectional  effect;  the  sulfur  should  reduce  jr-backbonding.  The 
C-0  bond  strengthens  as  the  chemisorptive  bond  weakens,  thereby  creating  a  high  CO 
frequency,  low  desorption  temperature  mode.  Although  two  such  modes  are  observed,  a 
third,  apparently  originating  from  increased  TT-backbonding,  does  as  well.  We  ambitiously 
hope  to  model  these  results  by  examining  the  effect  of  sulfur  coadsorption  on  the  binding 
of  CO  at  various  sites  and,  as  well,  gain  general  insight  into  the  mechanism  of  coadsorbate 
interactions. 

A  coverage  of  0=  j,  p(2x2)  per  adsorbate  species  was  chosen.  It  is  in  this  coverage  range 
that  the  three  C-O  stretch  peaks  first  appear  simultaneously  in  the  HREELS  specta®.  Each 
species  is  laid  out  in  a  square  net,  lattice  constant  a=4.98A.  The  calculated  band  structures 
of  CO22  and  S  unsupported  nets  axe  essentially  free  of  dispersion.  The  band  energies  are 
equivalent  to  the  molecular  orbital  energies,  thus  intraspecies  interactions  axe  negligible  at 
this  distance.  To  help  us  see  the  important  orbitals  of  both  adsorbates  and  their  relative 
energies,  we  show  in  Figure  1  the  total  density  of  states  (DOS)  of  the  bare  Ni  slab,  and 
alongside,  the  important  frontier  orbitals  of  S  and  CO. 


Figure  1  here 


One  of  the  prerequisites  for  the  electronegativity  arguments  concerning  poisoning  is 
that  the  electron  withdrawing  power  of  sulfur  be  observed  in  Ni(100)-p(2x2)S,  that  is, 
without  CO  coadsorption.  Ten  and  16  k  point  sets  were  used  for  substrate  plus  adsor¬ 
bate^)  systems  with  tetragonal  and  orthorhombic  symmetry  respectively  (see  Appendix 
for  calculations  testing  convergence).  The  sulfur  square  net  is  placed  on  the  substrate 
such  that  S  occupies  four-fold  hollow  sites  and  d(S-Ni)=2.19A30,  as  obtained  from  LEED. 
The  z  axis  is  consistently  oriented  perpendicular  to  the  Ni  plane,  and  the  x  and  y  axes 
parallel  to  the  surface  Ni-Ni  close  contacts.  The  charge  transfer  occurs  as  predicted;  the 
sulfur  adatoms  gain  0.756e~,  and  each  of  the  four  associated  surface  atoms  looses  0.355e~ 
relative  to  bare  Ni(100).  Although  the  charge  transfer  is  greatest  in  the  nickel  xz  and  yz 
levels  (— 0.115e~  each)  and  the  sulfur  z  (+0.481e~),  the  projected  DOS  indicate  that  the 
strongest  interaction  is  with  the  sulfur  x  and  y  rather  than  z  levels.  The  p  band  of  the  S 
square  net  lies  at  —  13.2eV,  well  below  both  the  d  block  and  the  Fermi  energy  of  the  bare 
substrate  (— 8.7eV).  The  large  dispersion  provided  by  the  interactions  with  surface  xz  and 


Energy  (eV) 


Figure  1.  In  the  left  panel,  the  bare  Ni(lOO)  DOS.  The  valence  orbitals  of  CO  and  S  are 
drawn  in  the  middle  and  right  panels  respectively. 
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yz  states,  such  as  the  representative  one  depicted  in  3,  forces  some  20%  of  the  S  x  and  y 
states  above  the  Fermi  energy  and  splits  the  major  S  x  and  y  peak  down  from  the  z  by 
~0.5eV.  Some  surface  xz  and  yz  character  appears  in  the  S  x  and  y  peak  near  —  14eV,  as 
seen  in  a  comparison  of  the  projected  DOS  integration  curves  in  Figure  2.  In  contrast,  the 
S  z  level  acts  essentially  as  an  inert  electron  sink. 

Figure  2  here 

The  total  DOS  (solid  line  in  Figure  2)  appears  very  much  as  a  simple  overlay  of  the  S 
square  net  DOS  (represented  schematically  on  the  right  by  bars  as  the  median  energies) 
and  the  bare  Ni(100)  DOS28.  The  binding  energies  of  the  sulfur  p  and  s  levels  agree  with 
those  obtained  by  UPS31  within  experimental  error.  The  S  adlayer  is  a  relatively  small 
pertubation  on  the  electronic  structure  of  the  Ni(lOO)  slab  as  a  whole,  but  does  cause 
substantial  charge  redistribution  and  rehybridization  at  a  local  level.  If  poisoning  does 
occur  via  the  substrate,  the  xz  and  yz  surface  levels  are  the  most  likely  mediators. 

We  consider  three  high  symmetry  sites  for  CO  coordination:  on-top,  bridging  and 
four-fold  hollow.  The  C-Ni  and  C-0  bond  lengths  axe  determined  by  LEED  as  1.80 A  and 
1.15A  respectively12*  for  on-top  coordination,  the  only  one  known  to  exist  independently.  In 
general,  our  computational  method  cannot  be  trusted  to  predict  the  energetics  of  adsorbate 
site  preference.  So  we  must  make  some  geometric  choices.  For  the  bridging  and  4-fold  sites, 
all  calculations  on  both  clean  and  sulfur  coadsorbed  surfaces  were  performed  at  two  limiting 
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Ni  surface  xz+yz 
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Figure  2.  The  solid  line  indicates  the  Ni(lOO)-p(2x2)S  total  DOS.  The  dotted  and  dashed 
lines  are  the  integrals  of  the  surface  xz  +  yz  and  S  x  +  y  projected  DOS  respectively. 
The  S  unsupported  square  net  DOS  is  represented  by  the  median  energy  bars  on  the 


geometries:  C-Ni  nearest  neighbor  distance  of  1.80A,  and  C-Ni  surface  distance  of  1.80A. 
The  trends  are  identical  but  consistently  more  pronounced  for  the  first  option,  which  we 
choose  to  present  for  the  sake  of  clarity. 

Characteristic  of  the  rather  weak  chemisorption  of  CO  onto  nickel  and  other  late  transi¬ 
tion  metals,  the  CO  molecular  levels  for  the  observed  on-top  geometry  remain  moderately 
discrete  and  energetically  stationary.  In  contrast,  one  third  of  the  bridging  and  4-fold  27t* 
levels  are  pulled  into  the  d  block.  Specific  features  of  the  DOS’s  will  be  discussed  later,  in 
conjunction  with  the  coadsorbed  systems. 

Our  models  adhere  closely  to  the  Blyholder  model  of  a-donation  and  7T-backdonation. 
The  electron  density  shifts  from  bare  Ni(100)  or  molecular  CO  to  Ni(100)-p(2x2)CO  are 
presented  in  Table  1.  The  extent  of  5 a  depopulation  falls  within  a  much  narrower  range 

Table  1  here 

than  does  the  amount  of  backdonation  into  27T*.  Bagus  et  al.S2a,i  feel  that  the  net  a 
interaction  may  contribute  little  to  the  chemisorptive  bond.  We  don’t  agree.  Rather, 
we  think  that  the  a  interaction  contributes  to  the  net  bonding  in  a  relatively  uniform 
manner  for  all  sites.  On  the  other  hand,  the  C-O  overlap  population  (o.p.)  follows  the 
2 ?r*  occupation  exactly  as  predicted.  Moving  through  the  coordination  series  of  molecular, 
on-top,  bridging  and  4-fold,  the  C-0  o.p.  falls  in  conjunction  with  growing  2 t*  occupation. 
The  CO  coordination  number  mediates  the  electronic  shifts  at  the  associated  surface  atoms 


Table  1.  Calculated  Results  for  Ni(100)-p(2x2)CO 


Atom(s)  or 

CO  Coordination  Mode 

Orbitals 

on-top 

4-fold 

bridging 

CO  total 

+0.249 

+0.777 

+0.638 

Electron 

CO  5a 

-0.379 

-0.368 

-0.398 

Density 

CO  2jr*(b> 

+0.372 

+0.694 

+0.598 

Changes^ 

Ni(, bridged 

-0.814 

-0.283 

-0.542 

Nii, unbridged 

-0.005 

— 

-0.002 

Ni„average 

-0.200 

-0.283 

-0.273 

Overlap 

c-ow 

1.044 

0.916 

0.942 

Populations 

C-Ni, 

0.847 

0.359 

0.631 

CO  binding  E  (eV)<d> 

-2.570 

-2.526 

-3.843 

•Relative  to  Ni(100)  or  molecular  CO.  Net  electron  gain  if  positive,  loss  if  negative. 


f 

' 

bOccupation  of  2jt*  and  other  degenerate  orbitals  are  given  for  each  individual  orbital. 
‘Molecular  C-0  o.p. =1.208 

dE(Ni(100)-p(2x2)CO)  -  E(Ni(100))  -  E(molecular  CO) 


the  largest  shift  coincides  with  the  lowest  coordination.  The  average  electron  loss,  however, 
correlates  well  with  the  gain  in  2ir*,  but  is  independent  of  5o.  These  two  effects  underline 
the  notion  that  the  chemisorptive  interaction  is  governed  by  the  extent  of  7r-backbonding 
against  a  fairly  uniform  a-donation  background.  Unfortunately,  we  cannot  compare  the 
strengths  of  the  chemisorptive  bonds  directly;  the  effect  of  coordination  number  on  overlap 
population  is  generally  nonlinear. 

Atomic  electron  transfers  must,  however,  be  used  with  care.  In  this  case,  both  CO 
interactions  work  to  decrease  the  surface  electron  density.  Most  of  the  d  block  lies  just 
below  the  Fermi  energy  so  that  some  states  will  be  pushed  above  it  by  dispersion  from 
any  source,  whether  from  below  (5o)  or  above  (2jt*)m.  The  two  effects  are  most  cleanly 
separated  at  the  on-top  geometry.  The  orbitals  able  to  interact  with  the  5 o  do  not  interact 
effectively  with  the  27r*,  and  vice  versa.  The  appropriate  surface  d  orbitals  (zJ  for  5a,  xz 
and  yz  for  2jt*)  experience  nearly  identical  charge  shifts.  Since  the  5a  occupation  is  very 
similar  at  the  three  sites,  we  can  safely  surmise  that  its  interaction  with  the  substate  is 
constant  and  the  total  surface  electron  loss  is  determined  by  the  2 jt\ 

With  regard  to  the  electronegativity  arguments  concerning  sulfur  poisoning,  we  suggest 
that  the  electron  withdrawing  power  of  CO  must  be  considered  as  well.  The  total  charge 
transfer  into  CO  is  approximately  one  third  of  that  into  the  S  of  Ni(100)-p(2x2)S  for  CO 
on-top,  but  nearly  equivalent  for  the  bridging  and  4-fold  coordinations.  Although  the  point 
must  be  reconsidered  for  the  coadsorbed  systems,  this  result  makes  rationalizations  based 


on  electronegativity  alone  less  self-evident. 


Four  Ni(100)-p(2x2)S+p(2x2)CO  systems  have  been  considered:  1)  CO  on-top,  4a,  2) 
CO  in  four-fold  hollow  retaining  tetragonal  symmetry  (4-fold  tetra),  4b,  3)  CO  in  four-fold 
hollow,  symmetry  reduced  to  orthorhombic  (4-fold  ortho),  4c,  and  4)  CO  bridging,  4d. 


4c  4d 


The  top  layer  of  the  three  layer  unit  cell  is  shown  in  the  schematics.  In  each  case,  we 
have  left  the  S  in  a  four-fold  hollow.  Behind  this  assumption  is  the  strong  propensity 
of  the  S  for  this  site;  to  our  knowledge,  no  deviation  from  such  adsorption  is  known 
experimentally.  CO,  on  the  other  hand,  is  known  to  have  low  barriers  to  moving  into 
alternative  coordination  modes. 

The  evolution  of  the  C-0  and  C-Ni,  o.p.’s  (Ni,=surface  atom)  upon  sulfur  coadsorp¬ 
tion,  shown  in  Table  2,  are  intriguing  in  light  of  the  Blyholder  and  “anti-Blyholder”  modes 


observed  by  Gland9.  Just  as  predicted  by  the  Blyholder  mechanism,  S  coadsorption  in 


Table  2  here 

the  4-fold  ortho,  4c,  and  bridging,  4d,  systems  induces  C-O  bond  strengthening,  C- 
Ni,  bond  weakening  and  reduces  the  CO  binding  energy.  Although  the  absolute  binding 
energies  calculated  by  the  extended  Huckel  method  are  unreliable  and  we  are  unable  to 
predict  adsorption  site  preferences,  the  relative  binding  energies  at  the  same  site  can  be 
used  with  more  confidence.  We  tentatively  associate  the  4-fold  ortho  and  bridging  sites  to 
the  high  frequency  C-O  stretch,  low  temperature  desorption  mode.  Exactly  the  opposite 
effects  are  observed  for  the  4-fold  tetra,  4b.  This  site  corresponds  to  the  “anti-Blyholder” 
low  frequency  C-O  stretch,  high  desorption  temperature  mode. 

Both  bonds  are  weakened  if  CO  is  bound  on-top,  and  the  CO  binding  energy  decrease 
so  sharply  that  the  Ni(100)-p(2x2)S+p(2x2)CO  system  is  less  stable  than  its  separated 
components.  The  anomalous  behavior  arises  directly  from  the  close  C-S  contact.  In 
this  geometry,  4a,  with  j  coverage  for  both  S  and  CO,  the  S-C(O)  seperation  is  only 
1.83A.  This  is  very  much  within  a  bonding  range;  indeed  typical  S-C  single  bonds  fall 
between  1.75  and  1.80A.  The  computed  C-S  o.p.  is  0.627,  which  is  the  largest  value  of 
any  site  (the  bridging  mode  is  a  distant  second  at  0.015).  The  surface  species  is  better 
described  as  molecular  SCO  than  as  coadsorbed  S  and  CO.  To  our  knowledge,  no  evidence 
of  SCO  formation  has  been  observed  for  any  coadsorbed  S/CO  surface.  The  energy  of 


Table  2.  Overlap  Populations  &  Binding  Energies 


CO  Site 

C-0  Overlap  Population 

C-Ni,  Overlap  Population 

ACO  binding 

Energy  (eV)^ 

Without  S 

With  S 

Without  S 

With  S 

on-top,  4a 

1.044 

0.929 

0.847 

0.579 

-3.391 

4-fold  tetra,  4  b  j 

0.916 

0.906 

0.359 

0.363 

+0.050 

4-fold  ortho,  4c 

0.916 

0.928 

0.359 

0.336 

-2.435 

bridging,  4d 

0.942 

0.974 

0.631 

0.581 

-1.256 

•(CO  binding  E  without  S)  —  (CO  binding  E  with  S) 


this  coadsorbate  geometry  is  very  high.  We  believe  it  is  energetically  inaccesible,  and  that 
what  we  are  in  fact  modeling  is  local  site  blockage. 

Essentially  all  other  calculations  to  date  have  been  performed  at  similarly  short  C-S 
contacts.  This  distance  is  crucial;  at  1.83A  we  observe  site  blockage,  at  2.66A  (4-fold  ortho) 
and  2.78A  (bridging)  the  Blyholder  result,  and  at  3.64A  (4-fold  tetra)  an  “anti-Blyholder” 
effect.  The  smallest  distance  must  be  classified  as  short-range  and  the  large  S-C  o.p. 
stongly  indicates  a  direct  interadsorbate  effect;  the  largest  is  clearly  long-range  and  its 
corresponding  o.p.  is  negligible  (—0.003),  pointing  to  a  surface-mediated  interaction.  The 
question,  then,  may  not  be  what  hi  the  absolute  range  of  interadsorbate  interactions  but 
rather  to  identify  the  type  of  interaction  active  in  each  range. 

An  important  point  in  our  analysis  is  that  the  2ir *  occupation  follow  the  strength  of  the 
chemisorptive  bond.  The  trends  shown  in  Table  3  are  consistent  with  the  x-backbonding 
mechanism;  invariably  C-Ni,  o.p.  reduction  is  coincident  with  2tt*  depopulation,  and  vice 
versa. 


The  relative  electron  withdrawing  power  of  the  two  adsorbates  is  reconsidered  in 
Table  4.  Relative  charge  transfer  is  defined  as  the  atomic  electron  density  in  Ni(100)- 
p(2x2)S+p(2x2)CO  minus  that  in  the  single  adsorbate  system  (Ni(100)-p(2x2)S  or  Ni(100)- 
p(2x2)CO).  A  positive  value  indicates  that  the  coadsorbed  species  aquires  electron  density, 
a  negative  value  indicates  a  net  loss.  According  to  this  indicator,  the  S  is  more  electroneg- 


Table  3.  Occupation  of  2 ir 


CO  Site 


on-top 
4-fold  tetra 
4-fold  ortho 
bridging 


Predicted  Changed 


Occupation 

Occupation 

Without  S 

With  S 

0.372 

0.515 

0.694 

0.717 

0.694 

0.686 

0.598 

0.521 

ative  than  4-fold  ortho  or  bridging  CO,  but  more  electropositive  than  4-fold  tetra  CO.  The 
results  are  consistent  with  the  notion  that  a  more  electronegative  species,  the  4-fold  CO 
tetra,  will  accept  more  backbonding. 


Table  4  here 

Although  the  results  thus  far  are  consistent  for  each  site,  the  question  remains  as  to 
why  these  sites  behave  in  fundamentally  different  ways.  Let  us  begin  with  the  bridging 
geometry.  The  arguments  for  the  4-fold  ortho  site  are  much  the  same. 

The  composite  characteristics  of  the  Ni{100)-p(2x2)S-|-p(2x2)CO  bridging  total  DOS 
can  be  seen  clearly  in  Figure  3.  The  DOS  of  the  CO  and  S  unsupported  nets  are  represented 
by  bars  at  the  right.  In  spite  of  energy  shifts,  their  peakedness  is  substantially  retained;  a 
minimum  of  50%  of  the  adsorbate  levels  fall  in  their  major  peaks  as  labeled.  As  previously 

Figure  3  here 

discussed,  the  direct  interaction  is  predictably  small  with  d(C-S)=2.78A  and  C-S  o.p. =0.015 
Some  mixing  occurs  between  the  CO  lx  and  the  sulfur  p  levels  since  they  are  energetically 
closely  matched.  But  the  lx  is  heavily  localized  at  the  oxygen  end  (see  Figure  1)  and  is 
therefore  relatively  inert  to  adsorption.  The  poisoning  must  occur  via  the  substrate. 

The  same  conclusions  can  be  drawn  from  the  electron  density  shifts,  of  which  a  selected 
few  are  listed  in  Table  5.  The  excess  charge  at  S  is  nearly  independent  of  CO  adsorption.  In 
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Table  4.  Electron  Transfer  to  Coadsorbates  Relative  to  Single  Species  Systems!11) 


CO  site 

S 

CO 

on-top,  4a 

-0.530 

-0.009 

4-fold  tetra,  4b 

+0.015 

+0.071 

4-fold  ortho,  4c 

-0.041 

-0.050 

bridging,  4d 

-0.051 

-0.179 

*A  positive  value  indicates  an  electron  density  gain  relative  to  single  species  adsorption, 
and  a  negative  value,  an  electron  loss. 
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CO  2tt* 


CO  5a 
S  3p 
CO  17T 


CO  4(7 


S  3s 


contrast,  the  local  surface  charge  is  very  much  determined  by  the  adjacent  surface  species. 
With  only  one  exception,  every  atom  of  the  initially  negatively  charged  surface  (bare 

Table  5  here 

Ni(100))  looses  electron  density  as  adsorbates  are  added  sequentially.  The  bridged  atoms 
of  the  coadsorbed  system  suffer  most  with  respect  to  either  single  species  slab. 

Let  us  focus  on  the  xz  orbitals,  which  lie  along  the  bridging  direction.  As  discussed 
earlier,  the  electron  withdrawing  power  of  S  is  felt  most  strongly  at  this  orbital  and  its 
degenerate  partner  in  Ni(100)-p(2x2)S.  The  largest  perturbation  of  CO  adsorption  onto 
either  the  clean  or  sulfided  surface  occurs  here  as  well.  The  xz  is  geometrically  well  suited 
for  interaction  with  both  the  5 a  and  the  2jt*x,  which  lies  along  the  bridging  direction, 
depending  on  the  phase  relationship  between  the  bridged  atoms.  The  out-of-phase  com¬ 
bination  can  participate  in  o-bonding,  5a,  and  the  in-phase  in  7r-backbonding,  5b.  The 
interactions  are  energetically  favored  as  well.  The  out-of-phase  states  will  predominate  in 


Table  5.  Electron  Density  Changes^,  Bridging  CO 


Ni(100)-p(2x2)CO 

Ni(100)-p(2x2)S+p(2x2)CO 

Atom(s)  or 

Orbitals 

relative  to 

relative  to 

Ni(100)  or  CO 

Ni(lOO)-p(2x2)CO 

Ni(l00)-p(2x2)S 

CO  total 

+0.638 

-0.179 

— 

CO  5 a 

-0.398 

-0.022 

— 

CO  27T’x(b> 

+0.746 

-0.144 

— 

Co  2jr*y 

+0.449 

-0.037 

— 

S  total 

— 

— 

-0.046 

S  x,y 

— 

— 

-0.028 

Surface  atom^ 

-0.542  (-0.002) 

-0.376  (-0.263) 

-0.563  (-0.090) 

s 

+0.004  (-0.005) 

-0.120  (-0.034) 

-0.039  (-0.001) 

z2 

-0.188  (+0.006) 

-0.024  (-0.065) 

-0.094  (+0.011) 

XZ<b> 

-0.262  (-0.087) 

-0.190  (-0.091) 

-0.337  (-0.063) 

yz 

+0.010  (-0.024) 

-0.109  (-0.089) 

+0.016  (+0.002) 

•Net  electron  gain  if  positive,  loss  if  negative. 
bBridging  direction  along  x. 

'Coordinated  (Uncoordinated) 
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the  lower  half  of  the  d  band  (near  5 a)  since  they  are  bonding  between  the  metal  centers. 
Likewise,  the  in-phase  set  will  be  concentrated  in  the  upper  half,  adjacent  to  the  2n* .  The 
interaction  with  the  2jt*x  is  particulary  strong;  the  occupations  of  the  originally  degenerate 
molecular  orbitals  differ  by  0.297e“  and  the  median  energies  are  split  by  ~leV.  Later,  we 
will  use  the  phase  relationships  between  the  Ni  atoms  and  a  technique  we  call  COOP  to 
identify  the  a  and  jt  interacting  states. 

The  stepwise  evolution  of  the  xz  can  be  traced  in  Figure  4.  The  integral  of  the  projected 
DOS  is  compared  for  Ni(100),  Ni(100)-p(2x2)CO  and  Ni(100)-p(2x2)S+p(2x2)CO,  and 
overlayed  onto  the  total  DOS  of  the  latter.  The  initial  localization  of  the  states  entirely 
within  the  d  block  is  disrupted  by  both  a  and  ir  effects.  The  mixing  is  evident  from 
the  substantial  redistribution  of  levels  into  both  the  5a  and  2n*  areas.  Both  effects  are 
magnified  with  addition  of  the  S  adlayer. 

Figure  4  here 

If  the  integrals  of  the  5 a  projected  DOS  are  compared  for  the  surfaces  with  and  without 
S,  we  can  see  that  sulfur  addition  pulls  some  of  the  5 a  states  into  the  CO  In-  and  sulfur 
p  peaks,  but  does  not  substantially  alter  either  the  median  of  energy  or  the  distribution 
of  levels  near  Cf.  On  the  other  hand,  a  similar  comparison  of  the  2;r*x  integrals  reveals 
that  the  27r*x  level  is  strongly  distorted  and  the  median  is  pushed  up  ~leV.  We  must 
remember,  however,  that  these  effects  cannot  be  caused  by  a  direc*  S-CO  interaction. 


yl 


Ni(100)-p(2x2)S+p(2x2)CO 
xz,  with  S  &  CO 
xz(  with  CO 
xz,  bare  surface 


CO  2n* 


Figure  4.  the  dashed,  dot-dashed  and  dotted  lines  are  respectively  the  integrals  of  the  xz 


projected  DOS  of  Ni(100),  Ni(100)-p(2x2)CO  and  Ni(100)-p(2x2)S-|-p(2x2)CO.  The 


solid  curve  is  the  total  DOS  of  the  latter.  The  median  energies  of  the  5 a  and  2tt* 


coadsorbed  levels  are  indicated  on  the  right. 


f  v.v. 


How  does  the  depopulation  of  the  CO  2tt*  come  about  on  the  sulfided  surface?  We 
saw  earlier  that  the  2 jt*  interacts  strongly  with  the  upper  part  of  the  metal  xz  band  (5b). 
The  S  p  orbitals  also  interact  with  the  same  part  of  the  band  and  push  it  up  (Figure  4). 
The  metal-2  jt*  interaction  becomes  stronger  since  the  energy  difference  is  reduced.  As  a 
result,  both  the  adsorbate  and  metal  levels  are  depopulated.  To  explain  why,  we  turn  to 
the  simple  two  level  interaction  diagram  used  so  commonly  for  molecular  systems,  and 
modify  it  to  indicate  the  inherent  level  dispersion  of  infinite  systems,  6.  The  schematic  is 


with  S 


no  S 


crude,  perhaps,  but  functional  nonetheless.  The  heightened  interaction  naturally  forces 
upward  the  resultant  antibonding  levels,  which  contain  both  metal  and  2x*  character. 
Contrary  to  molecular  systems  for  which  the  energy  of  the  HOMO  can  change  due  to 
interactions,  the  fermi  level  of  the  metal  surface  is  essentially  fixed,  regardless  of  the 
surface  species.  Our  three  layer  calculation  models  this  well;  the  greatest  shift  of  ef  after 


adsorption  is  0.13eV  from  the  value  of  the  bare  surface.  The  placement  of  ef  in  6  is  crucial. 


In  this  case,  the  arrangement  is  such  that  many  of  the  levels  are  pushed  above  ef.  A  similar 
interaction  in  molecular  systems  will  bring  about  a  gain  in  the  occupation  of  one  of  the 
levels  at  the  expense  of  the  other,  but  here,  because  ff  is  fixed,  both  metal  xz  and  CO  2n' 
levels  lose.  An  accompanying  diagram  could  be  drawn  for  the  CO  5 a  interaction  with  the 
bottom  of  the  xz  band  (5  b).  However,  most  of  the  levels  will  lie  much  below  so  that 
even  an  interaction  of  similar  magnitude  will  not  alter  the  5a  occupation  to  the  extent  of 
the  2jt*. 

Should  this  interaction,  by  which  S  modifies  the  ability  of  some  surface-localized  metal 
states  to  bond  to  the  2jt*  of  CO,  be  called  a  direct  one?  We  would  prefer  to  refer  to  it  as 
“through-bond  coupling” ,  a  concept  that  has  served  well  in  organic33®**  and  inorganic3*®*4 
chemistry.  Seemingly  localized  orbitals  (lone  pairs  in  two  parts  of  an  organic  molecule,  or 
two  metal  d  orbitals  in  a  binuclear  complex  connected  by  one  or  many-atom  bridges)  in 
two  or  more  parts  of  a  molecule  can  mix,  interact,  and  split  in  energy  as  a  result  of  mixing 
with  orbitals  of  intervening  atoms  or  groupings  of  atoms.  This  is  what  happens  here. 

% 

That  the  5a  interaction  is  modified  at  all  can  best  be  followed  indirectly  through  the 
Crystal  Orbital  Overlap  Population  (COOP)  curve35  of  the  bridged  nickel  atoms.  The 
COOP  is  created  by  weighing  the  DOS  by  its  contribution  to  the  Ni-Ni  o.p.  In  Figure  5, 
these  are  compared  for  the  clean  and  sulfided  systems.  From  5a  and  5b,  we  know  that 
some  Ni-Ni  bonding  states  (near  the  bottom  of  the  d  band)  can  interact  with  the  5a,  and 
that  some  antibonding  ones  interact  with  the  2n *.  Indeed,  the  resonances  with  the  CO 


are  picked  up  and  the  shift  from  bonding  to  antibonding  through  the  d  block  occurs  as 
expected.  However,  a  small  bonding  region  (circled  in  Figure  5)  appears  above  £f  of  the 
sulfided  slab.  Assuming  conservation  of  the  absolute  number  of  bonding  and  antibonding 
states,  this  packet  must  have  been  pushed  up  from  the  d  block  by  the  sulfur  and  again  by 
the  CO  5a.  Surface  levels  that  should  have  been  filled  by  o-donation  remain  empty.  The  a 
contribution  to  the  chemisorptive  bond  must  be  reduced,  regardless  of  the  stability  of  the 
5 a  occupation.  Since  the  number  of  affected  levels  is  not  substantial,  the  2jt*  occupation 
remains  the  overriding  factor. 

Figure  5  here 

Can  analogous  arguments  be  used  to  explain  the  anti-Blyholder  behavior  of  the  4-fold 
tetra  system?  The  charge  densities  computed  for  the  coadsorbed  slab  lie  much  closer  to 
the  single  species  values  than  do  their  bridging  counterparts.  A  few  are  listed  in  Table  6. 
That  the  magnitudes  are  small  should  not  concern  us  unduly  as  the  trends  tend  to  be 
more  reliable  in  similar  extended  Huckel  calculations.  Moreover,  the  trends  are  consistent 
among  themselves.  For  example  the  2 x*  occupation  remains  directly  proportional  to  the 
C-Ni.  o.p.  and  indirectly  proportional  to  the  C-0  o.p. 

Table  6  here 


The  4-fold  tetra  system  differs  from  the  others  in  several  respects.  The  coadsorbed 
4-fold  tetra  is  unique  in  that  S  and  CO  gain  electron  density  relative  to  the  single  species 
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Energy  (eV) 


Ni-Ni,  with  S  &  CO 
- Ni-Ni,  with  CO 


Figure  5.  The  bridged  Ni-Ni  COOP’a  for  the  Ni(100)-p(2x2)CO  and  Ni(100)-p(2x2)S+p(2x2)CO 

are  indicated  by  solid  and  dotted  lines  respectively. 
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Table  6.  Electron  Density  Changes^,  4-Fold  Tetra  CO 


Ni(100)-p(2x2)CO 

Ni(100)-p(2x2)S+p(2x2)CO 

Atom(s)  or 

relative  to 

relative  to 

Orbitals 

Ni(100)  or  CO 

Ni(100)-p(2x2)CO 

Ni(100) 

Ni(100)-p(2x2)S 

CO  total 

+0.777 

+0.064 

— 

— 

CO  5 a 

-0.368 

+0.009 

— 

— 

CO  2jt* 

+0.694 

+0.023 

— 

— 

S  total 

— 

— 

— 

+0.015 

S  x,y 

— 

— 

— 

Surface  atom 

-0.283 

-0.338 

-0.621 

-0.266 

2* 

-0.043 

-0.071 

-0.144 

-0.044 

xy 

-0.133 

-0.023 

-0.154 

-0.132 

xz,yz 

+0.003 

-0.129 

-0.126 

-0.011 

•Net  electron  gain  if  positive,  loss  if  negative. 


slabs.  It  is  also  the  only  system  for  which  the  modification  of  the  coadsorbates  on  the  Ni 
surface  is  nearly  additive.  The  sum  of  the  charge  transfer  between  the  bare  and  single 
adsorbate  Ni(100)  surfaces  is  very  close  to  the  shifts  between  the  the  bare  and  coadsorbed 
values.  The  largest  difference  is  in  the  degenerate  xz  and  yz  orbitals.  Individually,  S  and 
CO  remove  a  total  of  0.122e“  from  each,  but  if  coadsorbed,  0.014e~  more.  As  a  comparison, 
the  greatest  discrepancy  for  the  bridged  CO  systems  is  nearly  five  times  the  value  (in  z2). 
In  addition,  the  coadsorbates  work  on  entirely  different  parts  of  the  system.  This  can 
be  seen  by  comparing  the  relative  charge  transfers  between  coadsorbed  and  single  species 
slabs  to  coadsorbed  and  bare  slabs  (columns  2-5  in  Table  6).  The  effect  of  CO  on  the 
clean  or  sulfided  substrate  is  greatest  at  the  surface  xy,  and  smallest  at  the  xz  and  yz.  The 
reverse  is  observed  for  the  sulfur  modification  of  the  clean  or  CO  preadsorbed  slabs.  Thus 
the  adsorption  characteristics  of  either  adsorbate  are  nearly  independent  of  the  dispersion 
in  the  surface  caused  by  the  other.  The  same  information  can  be  obtained  by  overlaying 
the  appropriate  projected  DOS  curves  and  integrals. 

Clearly,  the  mechanisms  for  interadsorbate  interactions  are  severely  limited  in  this  sys¬ 
tem.  Not  only  does  the  distance  preclude  a  direct  effect,  but  adsorbate-surface  interactions 
are  segregated  such  that  the  potential  strength  of  an  indirect,  substrate  mediated  effect  is 
minimized  as  well.  At  this  point,  a  good  explanation  of  the  trends  is  still  obscure,  and  we 
hesitate  even  to  speculate  on  its  origins. 

The  effect  of  sulfur  on  CO  adsorption  appears  to  be  determined  by  the  interadsorbate 
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distance.  At  short  separations  the  S  and  CO  interaction  is  repulsive  —  we  see  site  blockage. 
At  intermediate  separations,  for  instance  in  the  energetically  favorable  geometry  where  S  is 
in  a  four-fold  hollow  and  CO  is  bridging,  we  see  the  workings  of  the  usual  Blyholder  model. 
The  2 jt*  of  CO  is  populated  less,  with  resultant  weakening  of  C-Ni  and  strengthening 
of  C-0  bonds.  We  trace  this  effect  to  the  modification  by  S  of  the  ability  of  specific 
surface  orbitals  (xs  and  yz)  to  backbond  to  2ir',  or  through-bond  coupling.  At  longer 
separations,  we  see  a  small  reverse  effect  (C-0  weaker,  C-Ni  stronger)  which  has  some 
experimental  support,  but  it  is  not  easy  to  rationalize.  In  general,  we  think  that  the 
simple  electronegativity  and  electron  transfer  explanation  of  coadsorption  effects  needs  to 
be  supplemented  by  a  detailed  orbital  analysis  in  terms  of  through-space  and  through-bond 
interactions. 
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Appendix 

The  calculations  were  performed  using  the  tight  binding  extended  Huckel  method210,6. 
The  H„’s  for  Ni  had  been  previously  determined28  by  charge  iteration  on  the  bulk  metal 
using  Gray’s  equations37.  The  Hu’s  for  C  and  O  had  also  been  determined  previously26  by 
three  cycle  iteration  on  CO  adsorbed  onto  a  four  layer  Fe(110)  slab.  A,  B,  and  C  iteration 
parameters  are  from  reference  38.  All  parameters,  including  the  uniterated  S  parameters, 
are  listed  in  Table  7. 

Table  7  here 

For  both  the  Ni(100)-p(2x2)S  and  c(2x2)S  systems,  the  Ni-S  bond  length  is  reported 
to  be  2.19A30.  The  Ni(100)-p(2x2)CO  geometry  is  d(Ni-C)=1.80A  and  d(C-0)=1.15A12c. 
Adsorbates  were  placed  on  one  side  of  the  three  layer  slab.  Ten  and  16  k  point  sets  were 
used  respectively  for  systems  of  tetragonal  and  orthorhombic  symmetry.  Convergence  was 
checked  for  10,  15  and  28  k  point  sets  on  the  Ni(100)-p(2x2)S  system.  The  maximum 
variations  among  these  sets  were  0.007eV  for  the  Fermi  energy,  O.OOleV  for  the  total 
energy,  0.010  for  orbital  electron  densities  and  0.001  for  overlap  populations. 


Table  7.  Extended  Huckel  Parameters 


Orbital  Ha(eV)  ft  ft  C2* 


Ni  4s 

-7.8 

2.1 

4p 

-3.7 

2.1 

3d 

-9.9 

5.75 

2.00 

0.5683  0.6292 

C  2s 

-18.2 

1.63 

2p 

-9.5 

1.63 

O  2s 

-29.6 

2.27 

2p 

-13.6 

2.27 

S  3s 

-20.0 

1.817 

3P 

-13.3 

1.817 

‘Contraction  coefficients  used  in  double  ?  expansion 
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Figure  Captions 


Figure  1.  In  the  left  panel,  the  bare  Ni(100)  DOS.  The  valence  orbitals  of  CO  and  S  are 
drawn  in  the  middle  and  right  panels  respectively. 

Figure  2.  The  solid  line  indicates  the  Ni(100)-p(2x2)S  total  DOS.  The  dotted  and  dashed 
lines  are  the  integrals  of  the  surface  xx  +  yz  and  S  x  +  y  projected  DOS  respectively. 

The  S  unsupported  square  net  DOS  is  represented  by  the  median  energy  bars  on  the 
right. 

Figure  3.  The  bridging  Ni(100)-p(2x2)S+p(2x2)CO  total  DOS;  major  peaks  are  labeled. 

The  bars  on  the  right  indicated  the  median  energies  of  the  CO  and  S  unsupported 
nets. 

Figure  4.  the  dashed,  dot-dashed  and  dotted  lines  are  respectively  the  integrals  of  the  xz 
projected  DOS  of  Ni(100),  Ni(100)-p(2x2)CO  and  Ni(100)-p(2x2)S+p(2x2)CO.  The 
solid  curve  is  the  total  DOS  of  the  latter.  The  median  energies  of  the  5a  and  2n' 
coadsorbed  levels  are  indicated  on  the  right. 

|  Figure  5.  The  bridged  Ni-Ni  COOP’s  for  the  Ni(100)-p(2x2)CO  and  Ni(100)-p(2x2)S-t-p(2x2)CO 

i 

i 

I  are  indicated  by  solid  and  dotted  lines  respectively. 
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